Granular instability in a vibrated U tube 
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We study experimentally the collective motion of grains inside a U shaped tube undergoing vertical 
oscillations, and we develop a very simple quantitative model that captures relevant features of the 
observed behaviour. The height difference between the granular columns grows with time when 
the system is shaken at sufficiently low frequencies. The system exhibits two types of growth: 
exponentially divergent (type I) and exponentially saturating (type II), depending on the size of 
the grains. The type I growth can be quenched by removing the air whereas the type II behavior 
can occur in the absence of air. There is a good agreement between the model proposed and our 
experimental results. 

PACS numbers: PACS 45.70-n, 45.70.Mg, 46.40.f 



Energy can be injected on a granular system using vi- 
brations to observe many interesting complcx phenomena 
like pattern formation [1|, segregation [3L internal collec- 
tive motion of the grains, among others [3J. Tubular con- 
tainers enhance the role of the walls in the interehange 
of energy and are extremely important from the indus- 
trial point of view. Many industrial processes involve the 
transport of grains through pipes and are unavoidably 
affected by mechanical vibrations, therefore their relia- 
bility can be greatly affected by the knowledge of the 
physics of vibrationally induced granular bulk flow Q . 

An interesting example of vibrationally induced granu- 
lar transport in a pipe can be realized in a partially filled 
U-tube. In this system vertical vibrations can cause the 
rise of one of the column of grains until the other one 
is emptied. The U-tube instability was first mentioned 
by Gutman in 1976 he pointed out that the inerease 
of the difference in level was due to the air pressure dif- 
ference aeting aeross the bottom of the vibrated tube. 
Ohtsuki et al. @ investigated the change in the level dif- 
ference for grains in two different size vessels conneeted 
at the bottom. They found that friction between parti- 
cles and the walls of the containers play a key role, and 
coneluded that for large particles the pressure difference 
due to the interstitial fluid is unimportant. 

More recently, other authors have considered the pres- 
sure difference due to the interstitial fluid as important 
for the appearance of the observed instability. Akiyama 
et al. 0] reported measurements of the pressure at the 
bottom of a partitioned bed linked near the base and 
vibrated at low frequencies. King et al. [İj], used a vi- 
brated U-tube filled with metal beads and water as an 
interstitial fluid and measured the pressure drop aeross 
the bottom of the container. The experiments reported 
so far are not conelusive, more research needs to be done 
to elucidate the mechanisms responsible for these insta- 
bilities and many open questions are being investigated 
[9-12]. Understanding the dynamics of vibrated granu- 
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FIG. 1: (a) Sketch of the experimental setup used to study 
the transport of grains in a vibrated U-tube; (b) Definition 
of the height difference Ah, and;(c) Movable counterweight 
used to align the system and maximize the verticality of the 
oscillations. 



lar bulk flow in these simple configurations may lead to 
improvements in the design of sand cores used in metal 
casting and the transport of granular materials through 
vertical pipes. 

Harmonic vibrations are typically used to fluidize a 
granular system and a fluidized state can generate flows 
or hydrostatic like forces. Reyes et al [13(, has shown 
that the observed fluid like behaviour can be in fact a 
very complex combination of jammed and air mediated 
fluidized states that appear during each eyele of oscilla- 
tion. This eyelic fluidization can produce interesting dy- 
namical features not yet understood, as for example the 
instabilities reported in this work. To explain our exper- 
imental results, we use the idea of eyelic fluidization in 
a very simple quantitative one dimensional mechanical 
model. We successfuly deseribe an exponentially diver- 
gent growth (type I) observed in a vibrated U-tube par- 
tially filled with grains. A similar model was previously 
developed to explain reverse buoyaney of a large sphere 
immersed in a vertically vibrated granular bed [14j . In 
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FIG. 2: (Color online) The curve on the right shows the ex- 
ponential growth of Ah with time, for glass spheres with 
d = (0.25 — 0.36) mm, vibrated at / = 8 Hz, and ampli- 
tude a y = 8.5 mm. The semilog plot in the inset shows the 
interval where the growth is exponential. Lissajous figures 
of the motion of the container corresponding to each curve, 
are shown on the right. The three curves to the left show a 
linear growth that results from the horizontal component of 
the driving oscillations. 



FigQ]we skecth the experimental setup, for the U tube 
experiments reported here. The tube has inner diameter 
D = 2 cm, separation between vertical arms of 6 cm and 
height of 14 cm. The vertical arms are linked at the top 
and connected to a vacuum system through two outlets. 
The tube is fbced to a plexiglas plate with a base attached 
to an electromagnetic shaker VTS-150. A function gen- 
erator provided a sinusoidal excitation signal a y sin(wi) 
with uj = 2ırf, where / is the frequency, a y is the vertical 
amplitude and t is time. 

To follow the height difference Ah defined in FigQJb) 
between the levels of the granular columns, we used a 
digital video camera and a stroboscopic light. To achieve 
maximum verticality in the vibration, we used a counter- 
weight at the back of the plexiglas base, in the form of a 
pendulum as shown in FigQJc). To monitör the vertical- 
ity of the vibrations we filmed the tube while a strobo- 
scopic light was flashed at a frequency slightly lower than 
the frequency of vibration of the shaker. Following a dot 
marked on the vessel we were able to obtain Lissajous 
figures for the motion of the container. In this way we 
could conveniently measure the horizontal amplitude a x 
and the phase shift <j) between the vertical and horizon- 
tal displacements. The grains used were spherical glass 
beads of different diameters d, between 0.18 mm and 0.36 
mm that we will refer to as small grains, with bulk static 
bed density p = (1.44 ± 0.04) g/cm 3 ; and poppy seeds 
of d w 0.8 mm that we will refer to as lar ge grains, with 
static bed density p = (0.70 ± 0.05) g/cm 3 . 

Type I growth was found only for small grains, vi- 
brated at low frequencies and specific values of ampli- 



tude so that ^ 0.02. In Fig|2] we show the type I 
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growth and the effect of the horizontal component in the 
behavior of Ah as a function of time, for glass beads with 
d = (0.25 - 0.36) mm, vibrated at / = 8 Hz. Different 
curves are produced changing the position of the mov- 
able counterweight shown in Figfljc) . Here the influencc 
of the horizontal component is clearly shown, the growth 
of Ah changes from exponential to a linear growth as 
a consequence of the horizontal component of the oscil- 
lations. The growth is exponentially divergent in the 
region between A and B. As the decreasing column be- 
comes empty, the curve deviates from exponential growth 
and eventually flattens (region B to C). 

As mentioned before, interstitial air is necessary for 
type I growth to occur. When the columns of small grains 
are set at different levels, and the air is evacuated, Ah 
decays exponentially to zero instead of growing (see Fig. 
[3]). Air, or more precisely, the insterstitial fluid, plays an 
important role in many other vertically vibrated granular 
phenomena [13-20]. 

When the U tube was partially filled with large grains 
the vertical vibrations produced a growth of Ah that was 
qualitatively different to the one observed for the small 
grains. An exponentially saturating growth, that we cali 
type II was obtained for these bigger seeds, as shown in 
Fig|4ja) . In this case the smaller column does not become 
empty. 

In order to further investigate the behaviour of the 
large grains in the absence of air and without an horizon- 
tal component in the vibration, we performed a molecular 
dynamics simulations of monodisperse and viscoelastic 
grains of diameter d p = 0.8 mm, with dissipative con- 
tact interactions [2Îj]. We used a system of N = 1700 
grains in a 2D U-tube of diameter d = 11.7 d p , vibrated 
at T = 1.64, without the presence of air and horizon- 
tal component. Results for the simulations are shown in 
Fig|4jb). For the lower curve the columns of the granular 
bed are initially leveled, so that the height difference be- 
tween the centers of mass of the two columns is Ah cm — 
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FIG. 3: The role of air is shown for glass spheres with d = 
(0.255 — 0.350) mm. (a) Exponentially divergent growth at 
low frequencies; (b) evacuation of the air up to a pressure of 
0.1 Torr, quenches the growth and the columns remain leveled 
when shaken at low frequency. If they are initially unleveled, 
Ah decays to zero as shown on the graph. 
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FIG. 4: (Color online) (a) Saturating growth for poppy seeds 
with d » 0.8 mm; (b) Simulations showing saturating be- 
haviour for particles with diameter d p vibrated at F — 1.64. 



d p . The low frequency vibrations produce results similar 
to those obtained experimcntally for poppy seeds. The 
top curve is produced by same granıüar bed, under the 
same conditions, except that initially the two arms are 
unbalanced so that Ah cm = 19 d p . Due to vibrations, 
Ah cm decreases to the same saturation vahıe of 11 d p . 
This implies that neither air nor an horizontal compo- 
nent is necessary for the saturating instability to occur. 

We believe that cyclic fiuidization is a key mechanisnı 
responsible for the growth of Ah in a vibrated U tube but 
there are two qualitatively different processes involved. 

A possible qualitative model for the saturating case can 
be described as follows: the vibrated granular medium 
acts like a granular bed that is thrown upward, detaches 
from the container and becomes fluidized, and during 
flight the bed feels a force due to an effective gravi- 
tational field pointing downward, acting as a restoring 
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FIG. 5: (Color online) Growth rate, for type I behaviour, as 
a function of the adimensional acceleration, for glass spheres 
with d = (0.18 — 0.21) mm. Three realizations at each value of 
F at a fixed frequency / = 12 Hz, changing the amplitud of the 
vibration. The solid line represents the model aproximation 
Eq. using the friction coemcient v as a free parameter. 



force. Later, during the settling time, in the same cy- 
cle of oscillation, the granular fluid changes to a jammed 
state by going through a compaction process that causes 
an increase in the interaction with the walls. This hap- 
pens in such a way that the largest column ends slightly 
higher than the shorter one, in each cycle. A quantita- 
tive model for this case will be presented in a future work 

The air mediated type I growth, reported here, corre- 
sponds to a process where the granular bed, significantly 
affected by the interstitial air, behaves similarly to a vis- 
cous liquid in a U tube, driven by an effective field pulling 
upward during part of the interval of time when the bed 
is in a fluidized state. The instability that leads to the 
type I growth can be described with a simple model sim- 
ilar to the one developed to describe the phenomenon of 
reverse buoyancy 14| . The U tube, partially filled with 
grains, has cross sectional area A = ırD 2 /4. The total 
mass of the granular bed is M g = pAL where p is the 
static bed density. The vertical position of the tube is 
given by w = a y cos(ujt). Whcn the system is fluidized 
we describe it as a liquid. We assume, for simplicity, 
that the granular bed fluidizes during an interval of time 
t that coincides with the portion of the period in which 
the acceleration given by the container is greater than 
gravity. This assumption could be relaxed to consider 
the actual time that the system fluidizes. If we define 
g e f = a y tü 2 cos(ujt) — g as the effective acceleration of 
the system, then we assume that the system is fluidized 
while g e f > 0. This effective acceleration will be approx- 
imately equal to the one felt by the granular columns, 
as long as the bed does not separates significantly from 
the bottom of the tube. The latter assumption holds for 
fine grains, at low frequency, in the presence of air. We 
take y in the growing side of the bed FigOJb) as the only 
degree of freedom necessary to describe the motion of the 
fluidized bed. In the reference frame of the container, the 
İD equation of motion for the evolution of y(t) is given 
by: 



M g y = m(y)g ef - vy, 



(D 



where m(y) — 2yAp is the mass of the unbalanced sec- 
tion (striped zone in FigfTJb)) and we have ineluded a 
dissipative term to account for friction losses. We inte- 
grate this equation över one period of oscillation to get 
the average displacement y, assuming, for this simplified 
case, that the integral is zero if evaluated outside the in- 
terval [-t/2, r/2], where t = - cos _1 (r _1 ). This gives 
us a simple differential equation for y. If we consider 
a vertically vibrated U tube with diameter D, partially 
filled with a granular medium of effective density p, solv- 
ing equation (1), we obtain an exponential dependence 
of y with time, for T > 1, with a growing rate 7 given by: 
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Using the friction coefficient v as the only free parame- 
ter we obtain a good agreement between experiment and 
model as can be seen in FigJŞ] The dependence on the 
adimcnsional acceleration F of the growing rate 7 given 
by equation (2) appears only within the square brackets, 
but p depends on the adimensional acceleration. This 
dependence can be incorporated by correcting for the ex- 
pansion of the granular bed as a function of T. Some 
complications may arise for some configurations because 
the expansion is not the same in the columns and in the 
link between them. The friction coefhcient may also de- 
pend on T. The dotted lines in Fig[5]were introduced to 
account for a 10% percent variation in the density that 
occur, when the granular bed is shaken 2l|. In our sim- 
plified model, we are not taking into account the fact 
that the actual fluidization time deviates from r with in- 
creasing T. At higher T the bed tends to remain fluidized 
even when g e f < 0, lowering y in each cycle, consequently 
lowering the value of 7. 

The model proposed is crude and could be improved by 
considering Gutman's model to incorporate the contribu- 
tion of the pressure difference across the bottom of the 
tube [5]. It seems plausible that the pressure difference 
across the bottom contributes to the transfer of grains 
between the two columns but it is not at ali clear how it 
can act as a force to drive the larger column up in such 
a way that the growth of Ah is exponentially divergent. 
Although the pressure difference across the bottom grows 
with the height of the rising granular bed, the mass of 
the same column also grows linearly with Ah, and the re- 
sultant forces oppose each other. In our model the force 
due to an effective gravity points upward and increases 
with increasing mass, producing an exponential growth. 

The fact that the range of parameters for which reverse 
buoyancy was reported [141 ] (fine grains, small frequencies 
and large amplitudes) is the same for which exponen- 
tially diverging growth was observed, together with the 
fact that the assumption of cyclic fluidization was suc- 
cessfully used to model both phenomena, gives support 
to the idea that cyclic fluidization is playing a key role 
in the dynamics of thesc interesting instabilities. 

We have reported an experimental study of a U-tube 
partially filled with grains undergoing vertical oscilla- 
tions. The evolution of the height difference between the 
levels of the granular material in each side of the tube 
shows distinct behaviors depending on grain size: type 
I, wich is an exponentially divergent growth and type II, 
wich corresponds to a saturating growth. Type II growth 
can occur in the absence of interstitial fluid, as was con- 
firmed by our computer simulations in a 2D U-tube. The 



type I growth observed experimentally was modeled as- 
suming cyclic fluidization, and relevant features of the 
observed behaviour were successfully captured by a very 
simple quantitative one dimensional model. 
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